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We demonstrate coherent tunable coupling between a superconducting phase qubit and a lumped el- 
ement resonator. The coupling strength is mediated by a flux-biased RF SQUID operated in the non- 
hysteretic regime. By tuning the applied flux bias to the RF SQUID we change the effective mutual 
inductance, and thus the coupling energy, between the phase qubit and resonator . We verify the modula- 
tion of coupling strength from to 100 MHz by observing modulation in the size of the splitting in the 
phase qubit' s spectroscopy, as well as coherently by observing modulation in the vacuum Rabi oscillation 
frequency when on resonance. The measured spectroscopic splittings and vacuum Rabi oscillations agree 
well with theoretical predictions. 



Superconducting qubit research has made tremendous 
strides in recent years. Superconducting qubits are rou- 
tinely made with coherence lifetimes approaching 1 /xs 
and beyond [IJ. Also, a number of coupled qubit exper- 
iments with fixed coupling between qubits have been per- 
formed [2-9|. Any real superconducting quantum com- 
puter, however, will be composed of an intricate network of 
many qubits coupled to each other in various ways, as well 
as coherent "quantum buses" that will manage the shut- 
tling of quantum information between distant qubits. This 
means that it will become increasingly difficult to imple- 
ment quantum information processing between many cou- 
pled quantum circuit elements with fixed coupling between 
elements. The need to control the coupling between var- 
ious elements, such as qubit-qubit interactions or qubit- 
quantum bus interactions is essential. A number of ways 
of implementing tunable coupling between quantum cir- 
cuit elements have been proposed in recent years [ 10l - IT4ll 
and performed experimentally |[T5l - [T9ll . One rather concep- 
tually simple way of implementing tunable coupling, pro- 
posed by [14J, involves use of a flux-biased RF-SQUID, 
operated in the non-hysteretic regime, as a tunable "flux- 
transformer" between elements. We have employed such 
a coupling scheme to coherently couple a superconducting 
phase qubit to a lumped element resonator. 

The circuit for this experiment is illustrated in Figure [T] 
It is composed of a phase qubit, with critical current /^o, 
shunt capacitance Cqs, and geometric inductance L^, cou- 
pled through a mutual inductance M^^, to the RF SQUID, 
referred to as "the coupler". The coupler has a critical 
current /co, geometric inductance Lc, and junction capaci- 
tance Cjc' It is coupled through a mutual inductance M^r, 
to the lumped element resonator of geometric inductance 
Lr and capacitance C^. All the junctions are via-style ion- 
mill junctions, and the capacitors were fabricated by use 
of "vacuum" capacitor technology [20|. There is also a 
residual mutual inductance Mq^ between the qubit and res- 
onator, which was gradiometrically designed to be as small 
as possible. 

The phase qubit is also coupled to external control and 



readout circuitry. A dc bias line, coupled to the qubit loop 
via a mutual inductance M^^, provides an external flux bias 
to the qubit. This bias controls the nonlinear Josephson in- 
ductance of the qubit that controls the energy level spacing 
between qubit states as well as level anharmonicity. The 
qubit is operated in a flux bias regime that creates an ap- 
proximately cubic metastable potential well of sufficient 
anharmonicity to reliably isolate the lowest two metastable 
states of the well L21jl. 

A microwave drive capacitively coupled via series ca- 
pacitance Cx provides the excitation energy to drive transi- 
tions between the two lowest qubit levels, labeled \g) and 
|e), respectively. A short 5 ns) measure pulse is then 
applied to induce tunneling of the |e) state to the adjacent 
stable well. The state of the qubit is read out via a DC 
SQUID coupled to the qubit's geometric inductance via a 
mutual inductance M^^q [22J. 
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FIG. 1 . (a) Circuit diagram for the phase qubit, coupler and res- 
onator. The qubit parameters are Iqo 0.6 /nA, Cqs — 0.6 pF, 
Lq 1000 pH, I3q 1.8, and Mqc ^ 60 pH. The coupler 
parameters are Ico — 0.9 /j^A, Lc ^ 200 pH, Cjc — 0.3 pF 
and /3c ^ 0.5. The resonator parameters are Lr ^ 1000 pH, 
Cr ^ 0.4 pF, and Mcr ^ 60 pH. (b) Optical micrograph of the 
circuit. 

The qubit's circulating current, Iq couples an amount of 
flux Mqclq into the coupler, generating a circulating cur- 
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rent Ic governed by the relation 



- sin {27T(j)x + 13 Jc 



(1) 



where ic = Ic/Ico is the normalized circulating coupler 
current, = + Mqclq) /$o is the net external flux 
applied to the coupler, and $c — 27rLc/co/$o < 1- This 
current then couples flux to the resonator via mutual in- 
ductance Mcr- For a given change in qubit flux, the flux 
change seen by the resonator depends on the offset current 
circulating in the coupler due to the external bias flux, 
The result is a tunable effective mutual inductance between 
the qubit and resonator given by 
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From equation [2] we see that the effective mutual induc- 
tance can be tuned anywhere between the following ex- 
trema: 
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corresponding to a null circulating current in the coupler. 
In particular, M^ff = 0, when the coupler circulating cur- 
rent is at the critical current. It is also worth noting that 
in the limit that /3c ^ 1, (Meff)max increases without 
bound. 

An interesting consequence of the changing effective 
mutual inductance between the qubit and resonator is that 
the resonator's frequency modulates with the applied flux 
as 
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where cj^q — l/\/IvCV. The measured resonator fre- 
quency is shown in Figure [2];b). 

We approximate the Hamiltonian of our system using the 
Jaynes-Cummings model in the rotating-wave approxima- 
tion, 



(5) 



where Hq — ^hujq&qz is the qubit Hamiltonian, Hr = 
[dral + 1) h(jJr is the resonator Hamiltonian, and the in- 
teraction term, Hj{^x) — T^9c{^x) {o'qCLr + o-~dl) de- 
scribes the exchange of a single excitation between the 
qubit and resonator at a rate proportional to 

CJr Mtotali^x) 
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where Mtotaii^x) = Meffi^^) + Mqr incorporates the 
direct mutual inductance between the qubit and resonator. 
The last two terms and H^, describe the coupling of 
the resonator and qubit to environments that give rise to 
the resonator decay rate and qubit decay rate 7 [23|. 

The lowest two levels of a the qubit and resonator form 
a four-dimensional joint Hilber space. We label the qubit's 
ground and first excited states as \g) and |e), respectively, 
and the resonator's ground and first excited states as |0) and 
respectively. According to equation [sj when the qubit 
is on resonance with the resonator, so that the detuning 
A = (jjq — (jJr = 0, individual eigenstates of the qubit and 
resonator, given by |^0), |eO), and |el) are no longer 
the eigenstates of the coupled system. The new eigenstates 
are found to be |^0), and |el) and the symmetric and anti- 
symmetric superpositions |±) = i |eO)) /\/2. The 
corresponding energy eigenvalues are Ego = O^E^i = 
2h(jJr and E± = h{ujr ± gd^x)). 

The coupler is first calibrated by sweeping its external 
flux bias, and measuring the effect on the tunneling 
probability of the \g) state of the qubit. By tracking the 
required applied qubit flux to maintain a constant total 
qubit flux (f)q = + MqJc) /$o such that the \g) state 
tunneling probability is approximately 10%, we can deter- 
mine the circulating current in the coupler as a function of 
Figure [2];a) shows the measured coupler circulating 
current as a function of applied coupler bias flux. 

The next step in the experiment is to demonstrate the ef- 
fect of the coupler on the quantum mechanical interactions 
between the qubit and cavity. We first look for a cavity in- 
teraction using well-established spectroscopic techniques 
[21 , 22 J. By use of figure [2];a) the coupler is set to the de- 
sired coupling strength and then qubit spectroscopic mea- 
surements are performed. When the qubit transition fre- 
quency nears the resonant frequency of the resonator, an 
avoided crossing occurs, splitting the resonant peak into 
two peaks. When the qubit frequency exactly matches the 
resonator's frequency (A = 0) the size of the spectroscopic 
splitting is minimized to g{^x) / ^- This whole cycle is re- 
peated for different flux biases applied to the coupler. We 
observe the size of the zero-detuning splitting modulate 
from a maximum of gmax/^ ^ 100 MHz down to no 
splitting (FigurejS] (a)). The spectroscopic measurements 
are a good indicator that the coupler is working, but we 
do not consider them to be proof of coherent coupling be- 
tween the qubit and resonator, because the length of the 
microwave pulse is longer (2^ 500 ns) than the lifetime of 
the qubit. 
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FIG. 2. Measurements of the dependence of Ic,(^r, and gc on ap- 
plied coupler flux, ^x/^o- The vertical dashed lines bracket the 
applied flux ranges for the waterfall data shown in Figur^ (a) 
The measured circulating coupler current as a function of applied 
coupler flux along with the theoretical fit giving /3c = 0.51. (b) 
Measured resonator frequency as a function of applied coupler 
flux, along with theoretical fit using /3c extracted from (a). The 
fit yields uJro/27T = 7.710 GHz. (c) Measured coupling strength 
as a function of applied coupler flux along with the theoretical fit 
using parameters extracted from the theory fits in (a) and (b). 

To demonstrate coherent tunable coupling we measure 
the vacuum Rabi oscillation period between the qubit and 
resonator as a function of applied coupler flux. Accord- 
ing to equation[5j neglecting dissipation, the zero-detuning 
probability of finding the system in the state |eO) as a func- 
tion of time is periodic and given by 

|(eO||*(t))|^ = |(e0|e^^^/^|e0)|2 

= ^(l + cos(2^,($,)t)). (7) 

When dissipation is added the oscillatory behavior of 
equation [7] decays exponentially with a rate given by 
lavg = (7 + ^)/2. When \4g{^x) \ = |^ - tI^ the system 
is critically damped and the |eO) state decays at the rate 
lavg with no oscillations [24J. The qubit's decay rate was 
measured to be 7 = l/Ti = 1/135 ns. The resonator's 
decay rate was not measured in this experiment directly, but 



previous experiments have found them to be smaller than 
K ^ 1/1000 ns [20|. If we take 7 >> a^, then the criti- 
cal coupling strength where the oscillations are expected to 
disappear is Qc/tv 1 MHz. 
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FIG. 3. Spectroscopic and time-domain data over the range 
^x/^o = -0.462 to ^x/^o = -0.366 bounded by the verti- 
cal dashed lines in Figure |2] (a) Waterfall plot of the spectro- 
scopic measurements of the |ib) states showing the splitting tran- 
sition from ^c(-0.462)/7r ^ 50 MHz through ^c(-0.421)/7r = 
to ^c(-0.366)/7r 40 MHz. The inset to the left is a 
3D plot of the qubit spectroscopy showing the avoided cross- 
ing transition through zero for applied coupler flux values close 
to ^x = —0.421. (b) The corresponding vacuum Rabi mea- 
surements demonstrating coherent modulation in the coupling 
strength ^c($a^). 

Experimentally, we excite the |eO) state by applying a 
short (Tp ^ 5 — 10 ns) pulse with the qubit on resonance 
with the resonator. The pulse is fast enough that the res- 
onator remains in its ground state during state preparation. 
We then measure the state of the qubit as a function of time. 
Figures [2jb,c) and [3] summarize the spectroscopic and 
time domain measurements. For g(^x) /tv > 10 MHz, 
the vacuum Rabi data are used to determine the coupling 
strength by applying a Fast Fourier Transform (FFT) to the 
measured probability data. For g(^x)/7T < 10 MHz, the 
FFT method is less reliable and the coupling strength is de- 
termined from the size of the splitting in the spectroscopy 
data at zero detuning. We clearly see the splitting in the 
spectroscopy shrink to zero when the coupler is "off", but 
the corresponding time-domain data do not appear to be 
exponential, as predicted by eqjsjwhen gd^x) — 0. There 
appears to be a rapid drop followed by a slow (2^ 7 MHz) 
oscillation in the data (Figure |4]). 

There is a higher order coupling channel not included 
in equation [6j resulting from the finite but small direct mu- 
tual inductance between the qubit and resonator, Mq^. This 
residual coupling strength is given by 123 
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For our design parameters Qresiduai ^ lOkHz, much 
too weak to account for the residual effect seen in the data. 
We believe the residual coupling effect is due to weakly 
coupled, spurious two-level system fluctuators (TLSs) in- 
teracting with the qubit at this frequency [21 J. We have 
used a scan of vacuum Rabi data that confirms these types 
of weak oscillations throughout the entire spectroscopic 
range, even at frequencies far detuned from the resonator. 
This indicates interactions with weakly coupled TLSs not 
seen in traditional spectroscopy measurements. Figure 
[4] compares the vacuum Rabi data taken at ^x/^o — 
—0.421 and the exponential and non-exponential Ti data 
taken at qubit frequencies far detuned from the resonator 
and where no TLS splittings were visible in the spectro- 
scopic data. 
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FIG. 4. A higher resolution trace of the occupation probability of 
the |eO) state when ^x/^o = -0.421 along with exponential Ti 
and non-exponential Ti measurements taken at a qubit frequen- 
cies largely detuned from the resonator. The non-exponential TI 
trace showed no evidence of a TLS interaction in the correspond- 
ing spectroscopy. 

We have demonstrated coherent tunable coupling be- 
tween a superconducting phase qubit and a lumped element 
resonator, using a separate, flux-biased RFSQUID as a me- 
diating element. Spectroscopically, the coupling strength 
was observed to modulate from a maximum 100 MHz 
to zero. The Vacuum Rabi oscillation frequency was ob- 
served to agree well with the spectroscopic measurements 
for ^ 7 MHz. The residual oscillations for 

weaker coupling strengths are believed to be due to spu- 
rious TLSs in the junction barrier and not the result of a 
residual coupling effect from the coupler. 
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